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Abstract 



Idealized merging argon plasma jets are simulated in 2D using both gas dy- 
namic and MHD models. Results indicate that peak pressures of several 
hundred kilobar can be achieved for high Mach number jets. Including a 
simple optically thin Brehmstrahlung radiation model and plasma targets 
shows that extremely high densities and magnetic fields can be achieved dur- 
ing jet merging on the order of 1000 times the initial density /field. Further 
^h investigations should include detailed ionization processes and more accurate 

radiation modeling to properly capture the radiation transport and subse- 
quent target compression. 

^H 

1. Introduction 
lO 

r** Recent interest in magneto inertial fusion [2J [51 E] has resulted in the in- 

vestigation of new schemes for magnetized target compression. The typical 

*— i foil liner approach suffers from the "stand off" problem where mechanical 

>■ connections to the device are destroyed with each shot. Many different solu- 

tions to the stand off problem have been described. One possible approach 
is the use of plasma jets for target compression J3J [21 HI El EJ- In this scenario 
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the target is moved into a chamber and a cylindrical array of plasma jets is 
fired simultaneously at the target. The jets merge to form a plasma liner that 
compresses the target. A new experiment at Los Alamos National Lab, the 
Plasma Liner Experiment (PLX) is currently under construction to explore 
plasma jet merging. This work is part of an effort to understand jet merging 
physics and develop numerical tools for modeling such a device. In particular 
this is a first effort at investigating jet merging with radiation. 

The current investigation uses ideal MHD for magnetized target simula- 
tions and the ideal Euler equations for simulations without magnetic fields. 
The MHD model used is single species so both the target and the jet are made 
of argon, in an ideal experiment the target plasma would be hydrogenic so 
a multi-species MHD model would be used for modeling. In the MHD case 
a simple optically thin Bremsstrahlung radiation model assuming constant 
ionization state is used. The effects of simple radiation model are so dra- 
matic that it becomes clear that precise knowledge of plasma chemistry and 
radiation transport will be very important in future theoretical/numerical 
investigations. Ultimately a collisional radiative model combined with a ra- 
diation diffusion is likely needed to make accurate predictions. 

The results presented show typical peak pressure and densities that could 
be achieved in a high powered merging jet experiment in the ideal case. In an 
actual experiment with similar jet parameters where jets are merging in 3D, 
peak pressures are likely to be significantly higher. The results also illustrate 
"radiative collapse" which could occur in this type of device if the jets were to 
remain optically thin during liner formation and collapse. The loss of thermal 
energy due to radiation in the jet during compression might mean that the 
ram efficiency of the jet remains high during the target compression phase 
resulting in efficient coupling of jet directed energy to target compression. 

2. Model 

Both gas dynamic (J2J) and magnetohydrodynamic (MHD) (J3J) models are 
used in this paper. The gas dynamic model is equivalent to the MHD model 
when the magnetic field is set to zero everywhere, however the gas dynamic 
model is used because it is faster computationally since there only 5 unknowns 
instead of 8. The gas dynamic problem is used in situations where there is no 
magnetic field. An ideal equation of state is assumed and a simple radiation 
model is used in several simulations. The radiation loss term is given in (fTl) 
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where the temperature T e is given in electron volts, n e is the electron number 
density in 7^/m 3 and Z e $ is the effective charge state. 
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Ideal MHD equations with a bremsstrahlung radiation loss term are used for 
problems with a magnetized target 
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The variable ^ is the correction potential. This term is used to propagate 
errors in the constraint equation V • B — out of the domain at the error 
propagation speed I\ This technique was originally developed in [T3]. The 
electron number density n e is computed from these equations as n e = Z eS — 
where m^ is the ion mass and the temperature is given by T e (kelvin) = 
k n(z +i) - ^ ne simulations performed in this paper are 2D meaning that 
the Z derivative is zero, however all terms including u z , B x and B y are kept 
even though their values remain zero throughout the simulation. Nautilus is 
a 3D code, however, 2D simulations were performed so that peak pressure 
convergence could be achieved. At low resolution the predicted peak pressure 
is dominated by numerical diffusion and is much lower than it should be. At 
a resolution of 800X800 cells or higher the numerical diffusion no longer 
dominates the measured peak pressure and convergence is achieved. This 
also suggests that in 3D, converged simulations would require 800X800X800 



/ 



cells and would take significantly longer to run. 

3. Numerical Approach 

Simulations were performed using Nautilus, a fluid plasma modeling code 
developed at Tech-X. In order to obtain jet propagation through vacuum a 
specialized approach was required. The approach taken is related to that 
used in (TTJ [12] and is described below. 

In these simulations the equations are solved at every time step using 
two different approaches. The first approach is a positivity conserving first 
order finite volume method. This technique maintains positive pressure and 
positive density in the solution even when there are density or pressure jumps 
of 9-12 orders of magnitude. The scheme is locally conservative and shock 
capturing, though only first order accurate. A first order positive solution 
throughout the entire domain is computed using this technique. This solution 
in cell i is called qf and the solution has associated with it one flux across 
each face s denoted ff. In this case s can be thought of as a global index of 
cell faces since each face is shared by two cells. 

To compute a higher order accurate solution a second scheme is used. 
The second scheme is a second order upwind MUSCL method which is also 
shock capturing and conservative, but not positivity conserving. The solution 
obtained using this approach in cell i is denoted qf and has associated with 
it one flux across each face s denoted /". 

Solving just using the first order accurate scheme results in solutions that 
are far too diffuse at the resolutions used in this paper and the jets expand 
much more rapidly than they should. As a result, it's necessary to use a 
combination of the first and second order scheme. 

Every face has associated with it one flux, the flux is added to one cell 
and subtracted from the other, this ensures conservation. The final solution 
is computed from a combination of fluxes from the positive solution and the 
accurate solution. The final flux f s chosen for face s becomes 
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# if P% > P and p% > p and P£ > P and p a L > p 



ff otherwise 



(4) 



where P R , p R , P£ and pi are the pressure and densities determined by the 
accurate scheme to the left and the right of face s. Pq and po are basement 
densities and pressures that are typically set to 1 x 1CT 6 to 1 x 10~ 9 times the 



peak initial pressure and density in the simulation. This approach ensures 
that the positive solution is used in cells where the pressure and density are 
predicted to be below the threshold. 

It's important to note, that if one simply sets a basement pressure and 
density (as is frequently done in MHD), mass, momentum and energy con- 
servation errors will occur that can drastically alter the results to these jet 
problems introducing unphysical effects that destroy the solution. This is es- 
pecially true when radiation losses are included. In our specific case, simply 
setting the background pressure and density to basement values introduced 
perturbations in the solution which caused small shock waves in our high 
speed jets and they lost symmetry before merging even occurred. This effect 
became exaggerated when the jets moved diagonal to the grid. Finally, if the 
background density and pressure were too high, say 1/1000 the jet density 
and pressure, the jets would create a large wake which would then interact 
with neighboring jets resulting in the development of shocks in the jet even 
before jet merging. 

4. Setup 

Parameter regimes for jet merging were explored by Cassibry[10] and by 
Parks [6J. Based on those results we've chosen jets that are in the fusion 
reactor regime. The PLX device is a spherical device where the jets are 
placed at a subset of the faces of a truncated icosahedron with up to 32 jets. 
In these 2D simulations we've used 6, 12 and 24 jets to gain some idea of the 
relationship between peak pressure and the number of jets in non-radiating 
simulations. In these gas dynamic simulation no target is included. The 
initial temperatures, jet diameters and pulse lengths are based on simulations 
suggested by CassibryfTU] for the PLX experiment. In the simulations that 
follow the following conditions are used. Jets are initialized to either Mach 
80 = 226 km/s or Mach 40 = 113 km/s. The 2D jet diameter is 8cm with 
a jet length of 50cm. The initial total temperature is 2eV while the initial 
number density is assumed to be 5.0 x lO 23 ^-. The plasma is an argon 
plasma and the bremsstrahlung radiation model assumes constant ionization 
state Z=9. Densities in the vacuum are set to 1 x 10~ 6 times the initial 
jet density and vacuum pressures are set to 1 x 10~ 8 times the initial jet 
pressure. The results have been split into two sets, the first set of results are 
gas-dynamic simulations where we are simply looking at the peak pressure 
and density ignoring radiation effects. These simulations help us bound the 



parameter regime for simulations with plasma chemistry. In the second set of 
simulations we use one jet configuration and look at the effect of our simple 
radiation model with a magnetized target. In the radiating case we've only 
included one jet configuration as that is sufficient to illustrate the important 
physical effect observed, namely the high peak density and magnetic field 
due to radiation losses in the merging jets and the target. Modifying the 
number of jets in the radiating case did not change the observed phenomena. 

5. Results 

Figure [T] shows results of a neutral fluid simulation at Mach 80 on a 
1600 x 1600 grid. Four snapshots are shown. The first sub plot (upper left) 
shows the jets just as they begin to merge and the liner is formed. Shocks 
develop where the jets begin to merge and are indicated by bright lines in the 
plot. The second sub plot (upper right) shows the jets after they've collided 
in the middle and the resulting shock wave begins expanding outwards. The 
third sub plot (lower left) shows additional shock expansion. Final sub plot 
(lower right) shows the development of Rayleigh Taylor instabilities. Figure 
[2] shows the same simulation using MHD with a 0.1 Tesla magnetized target. 
The first sub plot (upper left) shows the initial B z field. The second plot 
(upper right) shows the field near peak compression. The next two plots show 
the field again as it expands. Some Rayleigh Taylor instability is evident. 

Figure [3] and figure [4] show peak pressure, density and the total energy 
in the system for merging neutral jets at Mach 40 and 80 using 6, 12 and 24 
jets. No radiation was included in these simulations. Peak pressure reaches 
140kbar and density reaches 3kg/m 3 in the Mach 40 case and 800kbar and 
5.5kg/m 3 in the Mach 80 case. Total energy of the system is also plotted. 
Jets enter the domain through the boundary so total energy increases un- 
til the jets are turned off at which point the energy remains constant. As 
expected, pressure and density peak and then rapidly decline. Peak pres- 
sure is increased as the number of jets is increased in agreement with work 
performed by Cassibry [TO] . 

Figure [6] and [7] show plots of peak B z density, pressure and total system 
energy of Mach 40 and Mach 80 MHD jets with target a O.lTesla initial 
magnetic field. The initial target configuration is not in equilibrium and 
would expand if it were not for compression by the plasma jets. The field 
can be though of as an idealized solenoidal field. Results are shown with 
and without using the optically thin radiation model. Note that in [6] and 
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Figure 1: Simple jet merging using inviscid Navier Stokes equations. Snapshots of mass 
density. Simulation is 2D and we are looking down on plane of the jets. The first image 
shows liner formation as the jets begin to merge. Second image shows a series of Shockwaves 
as the jets interact at the center. Third image shows solution as the Shockwave expands 
outwards. 4th image shows the development of Rayleigh Taylor instabilities. 
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Figure 2: Simple jet merging using ideal MHD equations with a magnetized argon target 
(no radiation). Snapshots of B z . The first image shows the magnetic field in the target 
before the jets have begun to interact with the target field. The second image shows the 
target near peak compression from the jets. The third image shows expansion of the target 
and the fourth image shows the development of Rayleigh Taylor instabilities. 
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Figure 3: Merging Mach 40 gas dynamic jets without radiation. Total energy is recorded 
as MJ/meter since the simulation is 2D. Peak pressure reaches 140kbar in the 24 jet 
simulation. 



max pressure (kbar) vs time in n s 




Figure 4: Merging Mach 80 gas dynamic jets without radiation. Peak pressure reaches 
800kbar in the 24 jet simulation. 
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[7] the magnetic field and density have been plotted on a log plot. When 
radiation is included in the model the peak magnetic field reaches about 
100 times its value in the non-radiating case. What's more, peak field and 
density are maintained constant for several microseconds, corresponding to 
the length of time the jets are turned on. The results suggest that plasma jets 
might be able to be used used to obtain long lasting strong magnetic fields, 
furthermore they provide an upper limit on what to expect when complete 
radiation models are included (including reabsorption) . 

Using the Rosseland mean free path for bremsstrahlung radiation 

_ 64.0A 2 Tj /2 

L - 7 3 n 3 Cm ( 5 ) 



(defined in |7] for example) where A is the atomic mass number, T e is the 
electron temperature in keV and p is the mass density in -^-j, we can arrive 
at an estimate of the importance of absorption in different regimes. In these 
simulations peak plasma densities of around 5000 -§ were achieved while the 
initial jet was about .03— §. Figure 5 shows the Rosseland mean free path 
(in centimeters) for bremsstrahlung radiation as a function of density and 
temperature. When the mean free path is below about 1cm reabsorption 
will be significant. The plot suggests that reabsorption will be significant 
in the jet plasma during merging after it radiates a significant amount of 
energy away bringing it to around lOOeV in temperature. This plot assumes 
a constant Z e ff = 9. This also suggests that radiation from the magnetized 
target will be partially absorbed by the dense relatively cold liner. 
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Figure 5: Plot of bremsstrahlung Rosseland mean free path vs temperature and density 
assuming constant Z e tt = 9. In much of the plasma the optically thin model is accurate, 
however during liner formation the plasma will radiate significantly cooling it down at 
high density and decreasing the bremsstrahlung mean free path until a significant part of 
the radiation is re-absorbed. A radiation diffusion approximation can be used to capture 
this effect. This also indicates the importance of determining Z e ff accurately, as a larger 
Z e ff will decrease the mean free path and a smaller Z e tt will increase the mean free path 
by a significant amount. 
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Figure 6: Mach 40 MHD jet merging simulations with and without using optically thin 
bremsstrahlung radiation model. Log plots of the peak magnetic field show that the peak 
field reaches about 100 times the value achieved in the non-radiating case. 



6. Conclusion 

First steps in modeling plasma jet merging have been performed in 2D 
using ideal gas dynamic and ideal MHD equations with a simple optically 
thin radiation model assuming constant Z e ff. Gas-dynamic modeling with- 
out radiation or targets showed that merging jet peak pressures on the order 
of 100s of kBar could be achieved for the conditions simulated, peak pressure 
increasing with the number of jets. Results indicate late time Rayleigh Tay- 
lor instabilities in all cases and the heavy influence of radiation in the MHD 
simulations. Including radiation results in field enhancement 100 times that 
observed in simulations without radiation and resulting peak field compres- 
sion of several microseconds, substantially different than the non-radiating 
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Figure 7: Mach 80 MHD jet merging simulations with and without using optically thin 
bremsstrahlung radiation model. The field compresses to about 100 times the non- 
radiating case, similar to what occurs for the Mach 40 jets. 



14 



case. The radiation simulations illustrate the importance of properly cap- 
turing the argon plasma ionization state in addition to the importance of 
modeling the radiation properly. The optically thin model becomes invalid 
as the liner density rises and we would ultimately expect lower peak field 
compression as a result. The radiation in the jet effectively increases the jet 
Mach number, therefore increasing the ram efficiency of the jet. In these 
simulations, both jet and target were argon plasmas. In a real device, the 
target plasma would be deuterium, while the jets argon or some other high Z 
gas resulting in large radiation losses in the jet and lower radiation losses in 
the target. An accurate study will require complete atomic physics requiring 
multi-species MHD and radiation in both the optically thin and diffusion 
limit. These efforts will be part of future investigations. 
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